Description 

SYMBOLIC EVALUATION ENGINE FOR HIGH-PERFORMANCE 



SIMULATIONS 



Technical Field 

The invention relates to the field of computer simulation, and 
more specifically, a method of analyzing and predicting the behavior and 
performance of a system through a computer simulation. 

Background 

A simulation application provides a computer-based method of 
analyzing and predicting the behavior and performance of a system. A typical 
application will require a user to construct a model of the system by adding and 
connecting components. These components provide a mathematical 
representation of simple or complex physical phenomena in the form of 
equations. The application will collect and process the equations to a form that 
can be readily simulated using a numeric integration technique. 

Typically, time domain systems are used to model systems by 
performing a numeric analysis of the system, with very little manipulation 
performed symbolically prior to simulation. This set of equations is put together 
in the form of ordinary differential equations or differential algebraic equations. 
Unfortunately, because little manipulation is done symbolically prior to 
simulation, the simulation is not optimized for speed, hi other words, current 
systems fail to optimize for improved execution time, and in some cases may not 
be solvable at all. 

Constantinos C. Pantelides in "The Consistent Initialization of 
Differential-Algebraic Systems" describes a method of analyzing differential 
equations. The method identifies equations that constrain integrated variables. 



These equations are differentiated a number of times to determine the complete 
set of system constraints. While Pantelides simplifies the differential equations, 
the algorithm fails to identify variables that are to be constrained by the 
additional differentiated equations resulting from Pantelides. Sven E. Mattson 
and Gustaf Soderlind in "Index Reduction in Differential Algebraic Equations 
Using Dummy Derivatives" describe a method of dynamically selecting the 
constrained variables. The dynamic selection results in substantial performance 
degradation due to the overhead of the selection process and because the solution 
must be restarted whenever the constrained variables change. 

The prior art does not provide a simulation application that 
increases the ease of constructing a model, while minimizing the computational 
effort required to simulate it. The present invention is directed to overcoming 
one or more of the existing disadvantages in prior art systems. 

Summary of the Invention 

An embodiment of the present invention provides for a method of 
simulating a system. The method establishes equations modeling the system 
using terms having characteristics encapsulated within the term. Next, the 
method performs symbolic processing on the established equations for 
simplification. Additionally, the method performs system processing on the 
established equations for efficient simulation. 

Another embodiment of the present invention provides a machine- 
readable storage medium having stored thereon machine executable instructions, 
the execution of the instructions being adapted to implement a method of 
simulating a system. The method defines equations modeling the system using 
terms having characteristics encapsulated within the term. Next, the method 
performs symbolic processing on the established equations for simplification. 
Additionally, the method performs system processing on the established 
equations for efficient simulation. 



Additional embodiments of the present invention provide a 
method of symbolically processing a set of equations. The method includes 
assigning a portion of the set of equations to variables that have non-zero partial 
derivatives, and differentiating the remainder of the set of equations. Next, the 
method approximates an algebraic derivative for those equations that cannot by 
symbolically differentiated. The method includes symbolically integrating 
equations that cannot be assigned, followed by differentiating equations that add 
output derivates and integrating equations that add output integrals. Finally, the 
method eliminates an integral as each symbolically differentiated or integrated 
equation eliminates a numeric integration, such that the integral is converted to an 
algebraic variable by eliminating the derivative or integral relationship. 

Another embodiment of the present invention provides a machine- 
readable storage medium having stored thereon machine executable instructions, 
the execution of the instructions being adapted to implement a method of 
symbolically processing a set of equations. The method includes assigning a 
portion of the set of equations to variables that have non-zero partial derivatives, 
and differentiating the remainder of the set of equations. Next, the method 
approximates an algebraic derivative for those equations that cannot by 
symbolically differentiated. The method includes symbolically integrating 
equations that cannot be assigned, followed by differentiating equations that add 
output derivates and integrating equations that add output integrals. Finally, the 
method eliminates an integral as each symbolically differentiated or integrated 
equation eliminates a numeric integration, such that the integral is converted to an 
algebraic variable by eliminating the derivative or integral relationship. 

Additional embodiments of the present invention provide a 
method of eliminating an integral in a Pantelides algorithm. The method includes 
assigning a preferred integration location rank to one or more integrals. Utilizing 
the preferred integration location rank, integrals are assigned to equations. Next, 
the method eliminates the integration of assigned or solved integral variables. 
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[11] Another embodiment of the present invention provides a machine- 

readable storage medium having stored thereon machine executable instructions, 
the execution of the instructions being adapted to implement a method of 
eliminating an integral in a Pantelides algorithm. The method includes assigning 
a preferred integration location rank to one or more integrals. Utilizing the 
preferred integration location rank, integrals are assigned to equations. Next, the 
method eliminates the integration of assigned or solved integral variables. 

[12] Other embodiments of the present invention provide a method of 

block tearing equations. The method includes identifying variables in the 
equations in a block in which the variables appear linearly with constant 
coefficients. Next, the method solves nonlinear integration equations for their 
respective integrals, and determinines the solvability of the nonlinear equations. 
Next, the method solves the linear equations, and solves the nonlinear equations 
utilizing iterates and block variables solved from the linear equations. The 
method scans for solved for variables for identification of variables that are 
independent and may be removed from the block. 

[13] Another embodiment of the present invention provides a machine- 

readable storage medium having stored thereon machine executable instructions, 
the execution of the instructions being adapted to implement a method of block 
tearing equations. The method includes identifying variables in the equations in a 
block in which the variables appear linearly with constant coefficients. Next, the 
method solves nonlinear integration equations for their respective integrals, and 
determines the solvability of the nonlinear equations. Next, the method solves 
the linear equations, and solves the nonlinear equations utilizing iterates and 
block variables solved from the linear equations. The method scans for solved 
for variables for identification of variables that are independent and may be 
removed from the block. 

[ 1 4] Another embodiment of the present invention provides for a 

method of simulating a system. The method establishes equations modeling the 



system using terms having characteristics encapsulated within the term. Next, the 
method performs symbolic processing on the established equations for reducing 
the number of terms. Additionally, the method performs system processing on 
the established equations for efficient simulation. 

Another embodiment of the present invention provides for a 
method of simulating a component. The method establishes equations modeling 
the component using terms having characteristics encapsulated within the term. 
Next, the method performs symbolic processing on the established equations for 
simplification. Additionally, the method performs system processing on the 
established equations for efficient simulation. 

Brief Description of the Drawines 

The accompanying drawings, which are incorporated in and 
constitute a part of this specification, illustrate one embodiment of the invention 
and together with the description, serve to explain the principles of the invention. 

Figure 1 is a flowchart illustrating an overview of the simulation 
application for an exemplary embodiment of the present invention. 

Figure 2 is a flowchart of the symbolic engine for an exemplary 
embodiment of the present invention. 

Figure 3 is a diagrammatic illustration of a 2-D four-bar linkage 

model. 

Figure 4 is a flowchart of the redundant integral identification 
process for an exemplary embodiment of the present invention. 

Figure 5 is a flowchart of the system processing process of an 
exemplary embodiment of the present invention. 

Figure 6 is a flowchart of the block tearing process of an 
exemplary embodiment of the present invention. 

Figure 7 is a flowchart of the block sorting process of an 
exemplary embodiment of the present invention. 



Detailed Description 

Reference will now be made in detail to the present exemplary 
embodiments of the invention, examples of which are illustrated in the 
accompanying drawings. Wherever possible, the same reference numbers will be 
used throughout the drawings to refer to the same or like parts. 

An embodiment of the present invention provides for a method of 
simulating a system. The method establishes equations modeling the system 
using terms having characteristics encapsulated within the term. Next, the 
method performs symbolic processing on the established equations for 
simplification. Additionally, the method performs system processing on the 
established equations for efficient simulation. 

Figure 1 is a flowchart illustrating an overview of the simulation 
application for an exemplary embodiment of the present invention. The 
simulation 100 includes a number of processes that may operate sequentially or 
simultaneously depending upon their relationship as indicated in the figure. 
Detailed explanations of each process in the system will be provided at the 
appropriate portions of this specification. As an overview, at stage 105, 
equations are established using terms. Terms are newly defined building blocks 
for modeling a system. At stage 107, component equations are defined; at stage 
109 connectivity equations are defined; and at stage 111, user-defined boundary 
conditions are defined. At stage 113, symbolic differentiation and/or integration 
occurs, thus simplifying the equations for modeling during runtime. At stage 
1 15, user and component initial condition equations may be defined, establishing 
the initial condition of the system. At stage 1 17, numeric integration equations 
maybe defined and used in the establishment of the transient portion of the 
modeling of the system. 

Processing for the initial condition system and the transient system 
may be separated and performed in parallel processing, thus leading to 
efficiencies in modeling and execution. Thus stages 119, 121, and 123 of the 



initial condition system correspond in processing to stages 125, 127 and 129, 
respectively, of the transient system. At stages 1 19, 125 the initial condition and 
transient system receives their inputs. At stages 121, 127, system processing 
occurs. At stages 123, 129, the systems are ready for solution. Thus, systems are 
generated for simulation that are easier to construct and provide for a speedier 
and more easily repeated solution. 

Term building and equation definition will now be more fully 
described. A model is built by connecting components. A component 
mathematically represents some physical behavior by the equations associated 
with the component. The following example illustrates a component and it 
associated equations. For example, a component that models the translational 
acceleration of a point mass may be associated with the following equation: 

0=F-ma 

where F is the force applied, m is the mass, and a is the 
acceleration of the point mass. Here F and a are unknown variables and the 
causality of the equation has not been specified (i.e. is the equation solved for F 
or a). 

In the symbolic engine of an embodiment of the present invention, 
equations may be represented as a sum of terms: 

o = r,+r 2 +... 

A library of terms may be established, and may include basic 
terms, composite terms and collection terms. Basic terms include, for example, 
constants (e.g. 3.1) and scalars (e.g. 2x). Composite terms include, for example, 
exponents (e.g. T 2 ) and curve interpolation ( C(x,y) . Collection terms include, 
for example, summation (T } +T 2 +T 3 ) and products ( T, * T 2 * sin(T 3 ) ). Those 
skilled in the art will appreciate that other types of basic terms, composite terms, 
and collection terms could be used and are contemplated by embodiments of the 
invention. In addition, those skilled in the art will realize that other types of 
terms, other than basic, composite and collection, could also be used. 



Each term in the library may perform one or more mathematical 
operations on itself, such as, Evaluate, Simplify, Factorize, Differentiate, 
Integrate, and Invert. Other various mathematical functions can also be included 
in the library as the developer of a system may deem appropriate for 
implementation. 

Those skilled in the art will appreciate that this term approach 
simplifies the implementation of the symbolic engine as characteristics of a term 
class are encapsulated within the term and are transparent to the engine. The 
library can also be easily extended by defining new term classes to, for example, 
model the pressure drop across a hydraulic valve, or the body-to-global 
coordinate transformation matrix of a 3-D body. Because the library may be 
constructed by those having greater understanding of the physics of a system, less 
knowledgeable individuals may use such a library without the need to understand 
all of the underlying characteristics. In other words, those modeling a system 
need not have detailed knowledge of the operation of individual objects or 
components. 

In one embodiment, the terms are implemented in object oriented 
programming. In addition to the use of object oriented terms, embodiments of 
the present invention may also use a term group. Term groups enable multiple 
terms to be evaluated together. In some cases it may be more efficient or 
convenient to evaluate a group of terms together. For example, a component that 
interfaces a user-code function to produce three outputs (/, ,f 2 ,f 3 ) given four 
inputs (/, ,I 2 ,I 3 ,I 4 ), would be associated with the following term group of 
equations: 

0 = O 2 -/ 2 (/ 15 / 3J / 4 ) 
0 = O 3 -f 3 (I 2 ,I 3 ,I 4 ) 

Although the equations appear to contain three independent 
functions (/, ,f 2 ,f 3 ), they may be evaluated simultaneously as a term group 



when the user code is executed. The user code may be executed multiple times. 
However, for improved efficiency, the user code may not be executed multiple 
times. The symbolic engine handles such a case by using a new entity called a 
term group. A term group is a container that evaluates a specified function and 
stores the resulting values. These values are accessed in equations through an 
access term called a term group item. 

A flag may be set such that the function is evaluated only when 
the first term group item requests its value. Subsequent term group items access 
values stored in the term group following the initial evaluation. Using this 
approach prevents the function from being reevaluated for each item and allows 
the items to appear in the equations in any arbitrary order. For example: 

0 = O,- if(x > y) then T Gl (I, , I 2 ) else 2.0 

0 = O 3 ~T G3 (I 2 ,I 3 ,I 4 ) 

0^-7^,/3,/J 

In the above equations, upon a request to evaluate O x , if the 
condition in the first equation is false, the function in the term group need not be 
evaluated. However, when a request for evaluation of 0 3 is requested, T G3 
accesses its value, triggering the evaluation of all terms in the term group. On a 
later evaluation of the third equation, the value of T G2 may be recovered from 
storage within the term group. 

Figure 2 illustrates a flowchart of the symbolic engine for an 
exemplary embodiment of the present invention. At stage 205, the symbolic 
engine 113 assembles all the equations that collectively define the behavior of the 
model. These may include: the equations written by components contained in the 
model 107, connection equations that govern how components interact with one 
another 109, and boundary constraint equations specified by the user 111. The 
result may be a differential-algebraic system of equations (DAE) that can be 
simulated using numeric integration. The symbolic engine 1 13 performs 
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extensive processing on the DAE system to put it in a form that can be efficiently 
simulated. 

[43] The transient system of equations often has a DAE index (the 

number of times the equations maybe differentiated to allow for the solution of 
the derivatives) of greater than one. Directly solving such a system effectively 
determines some derivatives using numeric differentiation. This can result in 
significant loss of accuracy in simulation results or iterative convergence failure. 
Therefore, in one exemplary embodiment, the symbolic engine uses a modified 
form of the well-known algorithm by Pantelides to reduce the system to a DAE 
index of at most one. The algorithm is based on an efficient recursive scheme of 
assigning equations non-integrated variables. If an equation cannot be assigned, 
it is differentiated along with all the equations assigned to its incident variables. 
Three modifications have been made to the Pantelides algorithm, and are 
illustrated in Figure 2: 

[44] First, at stage 210, assignment of equations is restricted to 

variables that have non-zero partial derivatives. This improves the handling of 
discrete systems by limiting the propagation of the variable coloring scheme of 
Pantelides. For example, in the following set of equations: 
0 = 5,-^(0 

[45] 0 = S 2 - if (Si > 0) then 1 else 2 

0 = S 2 -f 3 (t) 

[46] The first and second equations are assigned S 1 and 

S 2 , respectively. The third equation cannot be assigned a variable and so it, and 
the second equation, are differentiated. Limiting assignment to variables with 
non-zero partial derivatives avoids the differentiation of the first equation. This 
is especially important if f x is not a differentiable function. 

t 47 ] Also, at stage 2 1 0, if an equation cannot be symbolically 

differentiated, an approximate/algebraic derivative equation is used. For example 
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if in the above set of equations f 3 is not differentiable, the derivative of the third 
equation can be approximated by: 

[48] 0 = S 2 -(S 2 -last(S 2 ))/h 

[49] Where h is the integration time step. 

[5 0] At stage 2 1 5, the above stages are repeated assigning equations 

non-differentiated variables. Equations that cannot be assigned are, if possible, 
symbolically integrated. This further reduces the number of variables that are 
numerically integrated. 

[5 1 ] At stages 220, equations are differentiated that add output 

derivatives. At stage 225, equations are integrated that add output integrals, hi 
addition to the equations differentiated and integrated in stages 210 and 215, the 
engine symbolically differentiates or integrates equations that provide more 
information to the user without increasing the number of numerically integrated 
variables in the model. These additional derivatives or integrals are referred to as 
output variables. For example, a simple gearbox writes the following equation: 

[52] Q = <d 1 -2co 2 

f 53 ] Where a> x and co 2 are the angular velocities of the two gears. If the 

derivative of a> x ( a, ) is defined in the model but not the derivative of co 2 (a 2 ), 
the above equation is differentiated to give: 

[54] 0 = a x -2a 2 

t 55 ] Here a 2 is added to the model, but is an algebraic variable that is 

determined by the above equation. Similarly, if the integral of co x {6 X ) is defined 
but not the integral of co 2 (9 2 ), the gearbox equation is symbolically integrated 
to give: 

[56] O = 0,-20 2 

t 57 ] At stage 230, as each symbolically differentiated/integrated 

equation eliminates a numeric integration, an integral should be eliminated, i.e., 
the integral may be converted to an algebraic variable by eliminating the 
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derivative-integral relationship. Pantelides does not identify which set of 
integrals should be eliminated. To determine the integrals to eliminate, the 
symbolic engine uses a new process that is based on the concept of a preferred 
integration location. A preferred integration location (PIL) is a hint from the user 
or component developer that given a choice, the integration of a specified 
variable is to be preferred over all other variables. 

Figure 3 illustrates a 2-D four-bar linkage model having an output 
arm 305, coupler arm 310, and crank arm 320. The model is located on a ground 
330. One of the four angles must be integrated and the remaining three are 
determined by algebraic constraints. Integrating the angle of the crank arm 320 
(0 4 ) is preferred as it will avoid singularities in the equations. Thus 6 4 is a 
preferred integration location. Each PIL has an associated rank. User specified 
PILs may have the highest rank. PILs specified by a component developer can 
have one of several lower ranks. All other potential integration locations have a 
default rank. The symbolic engine tries to assign to each integral equation (i.e. an 
equation that has been differentiated or is the result of symbolic integration) an 
integral that has the lowest PIL rank. The process is repeated allowing 
assignment to higher PIL ranks until all integral equations have been assigned. 
The integrals assigned to the integral equations are converted to algebraic 
variables. To avoid assigning linearly redundant variables, the linear terms in the 
integral equations maybe solved using Gauss Elimination. 

Figure 4 illustrates a flowchart of the preferred integration location 
process for an exemplary embodiment of the present invention. At stage 405, the 
process collects all integral equations and user and component defined PILs. At 
stage 410, the process identifies integral variables that appear linearly and 
nonlinearly in the integral equations. At stage 415, the current PIL rank is set to 
the default setting. At stage 420, each equation maybe assigned an integral that 
has a PIL rank less than or equal to the current value. The process may give 
assignment preference to integrals that appear linearly in the equation. 
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At stage 425, the process checks to see if all integral equations are 
assigned. If not, processing returns to stage 440 where the current PIL rank is 
increased. If so, processing continues to stage 430 where each equation is 
assigned an integral that appears linearly may be solved. The solved value is then 
substituted into other equations. If due to substitutions, an assigned variable is no 
longer in the equation, another integral, with minimum integration rank, is 
assigned to it. At stage 435, the assigned or solved integral variables are 
eliminated. 

As shown in process 115 of Figure 1, to start the simulation, the 
user may specify initial condition constraints. Most prior art simulation 
applications only allow users to set initial conditions on integrated variables. 
This limitation is caused by the DAE system being processed (for transient 
simulation) assuming that the integrated variables are determined by numeric 
integration. Thus, when solving for initial conditions, the only free variables that 
the user can constrain are the integrated variables. In some simulation 
applications, users may be able to set initial conditions on derivative or algebraic 
variables through the use of an indirect solution method, such as optimization. 
Such methods are computationally very expensive and prone to convergence 
difficulties. 

The symbolic engine of exemplary embodiments of the present 
invention allows users to constrain any system variable while directly solving for 
the initial condition solution. This flexibility is achieved by solving for the initial 
condition and transient solutions using two separate systems of equations that are 
independently processed. Both systems of equations contain component, 
connection, and boundary condition equations. Additionally, the initial condition 
system has user specified initial condition constraints, while the transient system 
has numeric integration equations. Another benefit of this approach is that it 
allows components to write different equations for the initial condition and 
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transient systems, e.g., a torque converter component may write an additional 
initial condition equation to specify the initial speed ratio across the converter. 
[63] Process 117 of Figure 1 handles numeric integration equations for 

the transient system. The symbolic engine of exemplary embodiments of the 
present invention may make no assumption about the numeric integration 
technique used for transient simulation. The integration technique adds equations 
to the transient system to constrain integrated variables. The type of equations 
added depend on whether the integration is explicit or implicit. An explicit 
integrator will add the following equation for each integral: 



t 64 ] 0 = x-k(h,x i _ 1 ,x i _ 2 ,...,x i _ 1 ,x i _ 2 ,...) 

t 65 ] An implicit integrator will add the following equation for each 
integral: 

[66] 0 = x ~ c ( h > x t-i > x i-2 >-,*,-, , Xt-2 ,».)- Hh,*i-1 , *,- 2 , x,_ 2 ,...)x 

t 67 ] where x is the integral, x is the derivative, and c and k are 



functions of the step-size (h) and the history of x and/or x. The symbolic engine 
processes the transient system equations with the added integration equations. 
For an implicit integrator, this inlines the numeric integration which results in a 
significant performance improvement. 
t 68 l Once the initial condition and transient system equations have 

been assembled, they are independently processed into a form that can be 
efficiently simulated. 

[ 69 1 Process 121 and 127 of Figure 1 handle system processing. Figure 

5 illustrates the system processing process of an exemplary embodiment of the 
present invention. At stage 505, the system equations are examined. At stage 
510, the system of equations may contain structural defects such as redundant 
equations or unconstrained variables. These defects are usually not unique; i.e., 
which equation is redundant or which variable is unconstrained depends on how 
they are selected. The symbolic engine identifies the defects such that users can 
most easily take action to correct them. The identification process is based on 
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assigning equations variables from a limited set. Multiple passes may be made 
increasing the set of variables that can be assigned, until no more assignments are 
made. The equations are sorted such that user-defined boundary and initial 
condition constraints are assigned last. Any remaining unassigned equations are 
redundant and unassigned variables are unconstrained. 

[70] The assignment passes made for the transient system may include: 

[71] 1 . Allow assignment to non-default boundary condition 

variables. Default boundary condition variables are specified by components. 

[72] 2. Allow assignment to any variable. 

[73] The assignment passes made for the initial condition system may 

include: 

[74] 1 . Allow assignment to non-integral and non-default 

boundary condition variables. 
[75] 2. Allow assignment to non-default boundary condition 

variables. 

[76] 3. Allow assignment to any variable. 

[77] Using this multi-pass approach, any unconstrained variables will 

be the ones users would expect to set boundary and initial condition constraints 
on. By sorting the equations such that user-defined constraints are assigned last, 
any redundant equations will most likely be identified as incorrect user 
constraints, which can be easily removed. 

[78] At stage 515, alias variables are removed. The engine identifies 

alias variables that need not be evaluated during the simulation. For example, in 
the following equation: 

[79] 0 = x + y 

[80] x can be considered to be an alias (scalar multiple) of y. By 

replacing all occurrences of x with -y, the evaluation and storage of x is avoided 
and potentially other system equations may be decoupled or further simplified. 
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At stage 520, the equations are partitioned. The symbolic engine 
transforms the system equations so that they can be efficiently solved during the 
simulation. The transformation rearranges the system into a series of subsystems 
or blocks that are to be solved in sequence. This process is commonly referred to 
as Block Lower Triangular (BLT) partitioning. The symbolic engine may use the 
efficient, well-known partitioning algorithm developed by Tarjan and 
implemented by Duff. 

Each block identified by partitioning is a coupled set of equations 
that may be solved simultaneously. For a lxl block i.e. a block consisting of one 
equation and one variable, the engine tries to symbolically solve the equation for 
the given variable. If the equation cannot be solved symbolically, a numeric root- 
finder such as Newton-Raphson is used to solve the equation during the 
simulation. For a NxN block i.e. a block consisting of N equations and N 
variables, a new tearing algorithm is used to determine how the block is solved. 
Tearing is a technique used to reduce the number of variables (iterates) that must 
be solved numerically. 

At stage 525, the block tearing process occurs. The tearing 
algorithm developed for the engine balances reducing the number of iterates with 
the ease of convergence of consecutive (at each simulation step) solutions. It 
assumes that in each iteration, the numeric root-finder solves a linear 
approximation of the system as represented by the equation residuals and 
Jacobian. The effectiveness of the each iteration is then a function of the range of 
the Jacobian i.e. the range of variable values (around the current value) for which 
the Jacobian is approximately constant. Hence, any equation that has constant 
partial derivatives should be solved symbolically as their range is infinite. 
Solving other equations could result in non-constant partial derivatives being 
multiplied, thus reducing the range over which the partial derivatives are 
approximately constant. 
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[84] Figure 6 illustrates the block tearing process of an exemplary 

embodiment of the present invention. At stage 605, the process identifies 
equations in the block in which the block variables appear linearly with constant 
coefficients. Terms in these equations that are not functions of the block 
variables may be functions of other variables in the system. These equations are 
referred to as linear equations, and all other equations are referred to as nonlinear 
equations. 

[85] At stage 610, the process solves nonlinear integration equations 

for their respective integrals. Solving these equations first reduces the coupling 
in the block equations. It is assumed that the highest order derivatives are 
iterates. 

[86] At stage 615, the process solves the linear equations using, for 

example, Gauss Elimination with pivoting. 

[87] At stage 620, the process determines the solvability of the 

nonlinear equations. A variable in a nonlinear equation can be solved for if: a) 
partial derivatives with respect to all incident variables are independent of the 
block variables; and b) the variable appears linearly in all nonlinear equations. 

[88] At stages 625, 630 and 640, the process solves the nonlinear 

equations given the current set of iterates and variables solved from the linear and 
integration equations. The process successively adds iterates (or tearing 
variables) at stage 640 until no more nonlinear equations can be solved. Tearing 
variables may be selected based on the following criterion (given by precedence): 

[89] - Cannot be solved for in any equation; 

[90] - Has continuous real values; 

[91] - Appears in the most number of equations i.e. couples 

equations; and 

[92] - Can be solved for in the least number of equations. 

[93] At stage 635, the process scans the solved-for variables to identify 

ones that need not be evaluated to solve the block. These independent variables 
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and their associated values are removed from the NxN block. They are added as 
individual 1x1 blocks to the sequence of blocks that make up the system. 

t 94 ] Returning to stage 530 of Figure 5, block sorting rearranges the 

system into a series of blocks. Partitioning guarantees that the blocks can be 
solved in the determined sequence, but this sequence is not unique. The symbolic 
engine takes partitioning a step further by sorting the blocks to minimize 
computational effort during simulation, while maintaining a valid evaluation 
sequence. The blocks are sorted into the following categories: static blocks, 
dynamic blocks, and output blocks. 

t 95 ] A static block evaluates a variable that has a constant value 

throughout the simulation. In an exemplary embodiment of the invention, these 
blocks may only be evaluated once, at the start of the simulation. They have the 
following characteristics: 

t 96 ] 1 • Symbolically solve for the block variable 

t 97 l 2. Are only a function of variables solved for in preceding 

static blocks 

[ 98 J Dynamic blocks have values that have an impact on numeric 

integration and dynamically change during the simulation. They solve for 
derivatives, integrals, and/or algebraic variables, symbolically or using a numeric 
solution technique. These blocks may be evaluated at each step during the 
simulation. They are sorted to move blocks down in the sequence until they are 
first needed. Consecutive independent blocks are rearranged to evaluate blocks 
with fewer iterates first. This avoids wasting effort on solving larger blocks 
(ones with more iterates), if the smaller block fails to converge. 

["] An output block evaluates a variable that does not impact the 

simulation or its results. These blocks are only evaluated when the user requests 
display of the respective variable value. Output blocks have the following 
characteristics: 

[ 1 00 J 1 • Symbolically solve for the block variable; 
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[101] 2. Block variable is not used by a dynamic block; 

[ 1 °2] 3 . Variable value can be evaluated post-simulation; 

[ 1 03 ] 4, Block variable cannot be required by a component during 

the simulation; and 

[104 J 5. Block variable cannot be an integral or derivative. 

[105] Figure 7 illustrates the block sorting process of an exemplary 

embodiment of the present invention. At stage 705, the process assumes that all 



blocks are dynamic. At stage 715, the block list is traversed from the top to the 
bottom to determine the static blocks, which are then removed from the list. At 
stage 720, the block list is traversed from bottom to top to identify output blocks, 
which are then removed from the list. At stage 725, the dynamic blocks are 
sorted. 

[106] Returning to Figure 5, stage 535 processes term compression. 

Term compression is the process of identifying duplicate terms to avoid their 
redundant evaluations. This produces a substantial simulation performance 
improvement by reducing the computational effort required to evaluate and solve 
the system equations. A compressor object implements the algorithm. The 
compressor contains a map of terms and their respective compression keys. A 
compression key is a string that uniquely identifies the value of the term. The 
compressor traverses from top to bottom the terms in the static and dynamic 
blocks. Each term, in sequence, is asked for its compression key. If the key is 
already in the compressor's map, the term is identical to a preceding term and so 
it is referenced to the preceding term's value. This avoids the evaluation of the 
duplicate term. If the key is not in the map, the key is added to the map if its 
evaluation is guaranteed. For example, in the following expression T x and T 2 are 
added to the map (if not already in the map), but not T 3 or T 4 as their evaluation 
is dependent on T 2 . 

[ 1 07] 0 = 7;- if(T 2 > 0) then T 3 else T 4 
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[ 1 08] Terms that are added to the map can also add auxiliary 

compression keys. The auxiliary keys are references to intermediate values held 
by the term. For example a curve interpolation term needs to determine the slope 
of the curve to calculate the term's current value. The curve term may hold on to 
this slope and add an auxiliary key for it. If subsequently the slope of the curve 
term appears in the system (such as in the partial derivatives required for solution 
of an iterate block), the auxiliary key is referenced to avoid re-evaluating the 
curve slope. 

Industrial Applicability 

[ 1 09] A hardware platform capable of implementing the system is now 

illustrated for one embodiment. By way of a non-limiting example, a system 
environment is described in which the features and principles of the present 
invention may be implemented. A system environment consistent with an 
embodiment of the present invention may include an input module, an output 
module, a computing platform, and a database. The computing platform may be 
adapted to include the necessary functionality and computing capabilities to 
implement simulations input through the input module and access, read and write 
to the database. The results of running the simulation are provided as output 
from the computing platform to the output module for printed display, viewing, 
or further communication to other system devices. Such output may include, for 
example, results of the simulation. Output from the computing platform can also 
be provided to the database, which maybe utilized as a persistent storage device 
for storing, for example, a library of terms. 

[110] In embodiments of the present invention, the computing platform may 

comprise a PC or mainframe computer for performing various functions and 
operations of the invention. The computing platform may be implemented, for 
example, by a general purpose computer selectively activated or reconfigured by 
a computer program stored in the computer, or may be a specially constructed 
computing platform for carrying-out the features and operations of the present 
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invention. The computing platfonn may also be implemented or provided with a 
wide variety of components or subsystems including, for example, one or more of 
the following: one or more central processing units, a co-processor, memory, 
registers, and other data processing devices and subsystems. The computing 
platform also communicates or transfers simulation output and input to and from 
the input module and the output module through the use of direct connections or 
communication links. 

[HI] Alternatively, communication between the computing platform and the 

input and output modules can be achieved through the use of a network 
architecture. In an alternative embodiment, the network architecture may 
include, alone or in any suitable combination, a telephone-based network (such as 
a PBX or POTS), a local area network (LAN), a wide area network (WAN), a 
dedicated intranet, and/or the Internet. Further, the network architecture may 
include any suitable combination of wired and/or wireless components and 
systems. By using dedicated communication links or a shared network 
architecture, the computing platform may be located in the same location or at a 
geographically distant location from input module and/or output module. 

[112] The input module may be implemented with a wide variety of devices to 

receive and/or provide the data as input to the computing platform. The input 
module may include an input device, a storage device, and/or a network. The 
input device may include a keyboard, a mouse, a disk drive, video camera, 
magnetic card reader, or any other suitable input device for providing system 
information or modeling information to the computing platform. Memory device 
may be implemented with various forms of memory or storage devices, such as 
read-only memory (ROM) devices and random access memory (RAM) devices. 
The storage device may include a memory tape or disk drive for reading and 
providing, for example, simulation data or term library information, on a storage 
tape or disk, as input to the computing platform. 
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[113] The output module may include a display, a printer device, and/or a 

network interface for receiving the results provided as output from the computing 
platform. As indicated above, the output from the computing platform may 
include one or more results of simulations. The output from the computing 
platform may be displayed or viewed through a display (such as a CRT or LCD) 
and printer device. If needed, a network interface may also be provided to 
facilitate the communication of the results from the computer platform over a 
network (such as a LAN, WAN, intranet or the Internet) to remote or distant 
locations for further analysis or viewing. 

[ 1 1 4 J The preceding text describes a comprehensive symbolic engine 

that assembles and processes a system of equations into a form suitable for 
numeric simulation. The engine implements several new processes to maintain 
ease of model construction while enhancing simulation performance. 

t 1 1 5 ] It will be readily apparent to those skilled in this art that various 

changes and modifications of an obvious nature may be made, and all such 
changes and modifications are considered to fall within the scope of the appended 
claims. Other embodiments of the invention will be apparent to those skilled in 
the art from consideration of the specification and practice of the invention 
disclosed herein. It is intended that the specification and examples be considered 
as exemplary only, with a true scope and spirit of the invention being indicated 
by the following claims and their equivalents. 



